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Chapter 1 
INTRODUCTION 

1.1 BACKGROUND 

In recent years arcjet thrusters have found increased use in aerospace applications. 
Originally developed in the 1950s and 1960's, the arcjet thruster can be described as a 
device in which an arc is created in the presence of a propellant, with the resultant hot 
gases expanded through a nozzle creating a thrust in the opposite direction. The cutaway 
view of an arcjet thruster is shown in Fig. 1.1. The thruster has been described in detail 
elsewhere [1]. 


INJECTION 



Fig. 1.1 Cutaway view of an arcjet thruster (Source: [4]) 
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Because of lack of sufficient space electric power (such as solar panels) and 
forseeable near term applications, work on the thruster was temporarily abandoned in the 
1960’s. However, new space mission requirements along with major advances in power 
electronics technology over the last two decades have caused a dramatic revival of interest 
in the arcjet thruster [2]. The primary interest in arcjet thrusters stems from the fact that 
they are up to two times more efficient than present chemical engines and resistojet 
technologies [3]. The resultant savings can be used to realize: 1) increased satellite 
lifetime, 2) larger satellite payloads, 3) use of a smaller rocket booster, or any 
combination of these [3]. Currently, thrusters of various power levels are being 

researched primarily for prime propulsion and North-South/East-West stationkeeping of 
geosynchronous satellites [1, 2, 4, & 5]. 

The standard output of a satellite solar panel is approximately 30 volts DC, while 
the arcjet thruster requires an operating voltage greater than or equal to 100 volts DC, and 
an operating current of approximately 10 amps DC. It is therefore obvious that a DC to 
DC converter will be required by this system. The convener will also serve to limit the 
following: 1) arcjet voltage fluctuations caused by variations in the solar panel output 
voltage (voltage regulation); 2) solar panel output voltage fluctuations caused by 
variations in the arcjet voltage (input/output isolation). Further details concerning 
converter design can be found elsewhere [2, 4]. 

It has been shown that the arcjet thruster has a static volt-ampere (I-V) 
characteristic which exhibits a decrease in voltage with increasing current (negative 
resistance slope) over the practical operating range [4]. The static volt-ampere 
characteristics of an arcjet thruster for various gas flow rates is shown in Fig. 1 .2. 

Once the required converter has been assembled, it must be tested to see if it will 
perform as desired. By connecting resistors to the converter output, resistive load test can 
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Mass Row kg/sac 
O 4.8x10 * 5 
□ 4.5x10 ' 5 
A 4.3x10 * 5 


Fig. 1.2 Static Volt/ Ampere characteristics for various flow rates (Source: [4]) 


be performed to determine load and line current regulation, load transient response, output 
ripple and power efficiency [2, 5]. The resistive load tests provide helpful information 
needed to evaluate converter performance. However, they provide no indication as to 
whether the converter will be stable when connected to a load device which has the 
negative resistance I-V characteristic shown in Fig. 1.2. 

One of two methods can be used to indicate whether the converter will be stable 
when operating with a negative resistance load. The converter can be tested either by 
using the actual load device, or by using a device which simulates the I-V characteristics 
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of the actual load 

Using the actual load device for converter testing has several drawbacks. First of 
all, the converter designer may not have access to the necessary equipment required for 
arcjet thruster operation, such as a vacuum chamber, fuel (or propellant), or even the 
thruster itself. Acquisition of the necessary fixtures and equipment would be time 
consuming as well as expensive. It might be possible to test the converter at a facility 
which has access to the required equipment, but this to can be costly and time consuming 
as, the nearest facility could be hundreds of miles away. Test complexity can also 
present a problem as the designer must become familiar with the necessary procedures 

required when operating a thruster, such as vacuum chamber operation, proper fuel 

\ 

mixtures and flow rates, and arc establishment criteria. 

By using an electronic load which simulates the I-V characteristics of an arcjet 
thruster, all of the problems mentioned above can be eliminated 


1,2 OBJECTIVE 

This study considers the development of an electronic load (hereafter referred to 
as a simulator) which has an I-V load characteristic similar to that of the arcjet thruster, 
shown in Fig. 1.2. The simulator was constructed in the laboratory using relatively 
inexpensive power electronics components. To avoid unnecessary labwork in future 
applications, care was taken so that the simulator adequately resembled the load and was 
free of internal instabilities. After construction and successful operation, the simulator 
was then used to provide an experimental demonstration of the stability of various arcjet 
supply arrangements. 
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To evaluate simulator performance, the simulator was tested with an actual current 
mode controlled (CMC) converter designed for a 1 kilowatt arcjet thruster. Upon 
successful operation with the simulator, the converter was then tested on an actual arcjet 
thruster at the NASA Lewis Research Center in Cleveland, Ohio, after which, the static 
and dynamic load responses of the simulator and thruster testing were compared. 

In addition, methods of varying the simulator dynamic impedance so that it could 
be made to more closely resemble that of an arcjet thruster, were also explored. Finally, 
possible simulator applications with other arc producing devices (such as in welding) 
were investigated. 



Chapter 2 


SIMULATOR DEVELOPMENT 

2.1 CIRCUIT DEVELOPMENT 

2.1.1 Introduction 

Upon further examination of Fig. 1.2, it can be seen that the static I-V load 
characteristic of an arcjet thruster is similar to that of a constant-power (IV-k, where k is 
a constant) curve. The dynamic impedance data from [4] also suggest that the arcjet has a 
load characteristic similar to IV-k for frequencies up to about 60 kHz. Therefore, if the 
simulator is constructed in such a way that it produces a constant-power I-V load 
characteristic, it will then have an I-V load chareteristic similar to that of an arcjet 
thruster. This leads to the development of the two-loop feedback system shown in Fig. 
2 . 1 . 

Examination of Fig. 1.2 reveals that the large signal, steady state operation of an 
arcjet can be described by the function: 

v«/*( 0 

Since for large signal, steady state operation the compensation capacitor (C fl ) of 
Fig. 2.1 will block current flow through its path, the multiplier output voltage (v a ) must 
equal the reference voltage (v re j). Thus, 

V ref =V a 

where for the particular multiplier used (MC1494L) it can be shown that v a of Fig. 2.1 is 
given by: 

v a = k i + V + *3' v + i AV‘ ,v = V 
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Fig. 2.1 Initial cicuit design 

Rearrangement of (2.1a) yields: 

V = Where ^ = k i‘ k v‘ k M ( 2Ab) 

Originally £, and kj were set to zero, such that the large signal, steady state operation 
of Fig. 2.1 would be determined from, 

(Z2a> 

or, 

i*v = k where k = - 7 ^ ( 2 . 2 b) 

k r 

as desired. 

The circuit in Fig. 2.1 can be subdivided into the following three parts: 1) the 
power multiplier, 2) the PI compensator, and 3) the final stage. 

2.1.2 The Power Multiplier 

The purpose of the power multiplier is to produce a constant power I-V 
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characteristic at the terminals of the final stage by multiplying a sense of the final stage 
voltage by a sense of the final stage current, and providing this signal as an input to the PI 
compensator. 

2.1.3 The PI Compensator 

The purpose of the PI compensator is to provide enough gain such that sufficient 
simulator bandwidth is achieved, while at the same time providing a means of 
compensation to aid in stabilizing this two-loop feedback system. 

/ 

2.1.4 The Final Stage 1 

The desired I-V characteristics should be produced at the terminals of the final 
stage when they are connected to an appropriate arcjet power supply. Almost all of the 
power dissipated by the simulator (1 kW for this design) will be dissipated by the final 
stage. Therefore, it will be necessary to parallel several transistors when constructing the 
final stage transistor (Ql) of Fig. 2.1, in order to keep from exceeding specified device 
ratings on current, power, and operating temperature. The emitter resistance (R E ) is used 
to ensure that the final stage transistor current divides equally among each of the 
paralleled transistors. 

The dissipative resistance (R D ) is used to reduce the amount of power dissipation 
in the final stage transistor by an amount equal to i 2 -R D , thereby allowing the use of fewer 
transistors and a smaller heat sink. For this design i will be approximately 10 A. 
Therefore, a relatively low value of R D will provide sufficient power reduction in the final 
stage transistor. It should be noted that the final stage can be incrementally approximated 
by a current-dependent voltage source in series with R D (see Fig. 2.2). From simple 
voltage division it can be seen that: 
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where Rsource represents the internal impedance of the arcjet power supply connected at 
the final stage terminals. Therefore, as long as R D « Rsource then V terminal ~ f° r all 

Rsource • Thus , the low value of R D also aids in producing a simulator terminal response 
that is independent of the impedance of the arcjet power supply. 


The final stage transistor was constructed from Insulated Gate Bipolar Transistors 
(IGBT’s). Reasons for choosing IGBT’s include simplicity of gate drive, low thermal 
resistance per unit cost, and high power handling capabilities. Gate resistors (R c ) were 
connected to the gates of each IGBT to prevent parasitic oscillations in the paralleled 
transistors, as described in [6]. 


As will be seen shortly, the feedback resistance (Rj) transforms the final stage into 
a transresistance amplifier, converting it to the current-dependent voltage source 
mentioned above. The resistance R x was added to form the current drive for the 


transresistance amplifier. 
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2.2 THE FINAL CIRCUIT 

The final working version of the simulator is shown in Fig. 2.3. Each R E was 

chosen to be 4.7 Q to ensure equal current division among the 16 IGBT’s and to provide a 
sufficient current sense signal to the current sense amplifier (OP-27). The input voltage 
range of the multiplier (MC1494L) is ±10 V. As a result the voltage divider and the 
OP-27 are arranged such that the simulator output voltage and current can go as high as 
200 V and 20 A, respectively, before the multiplier inputs are overloaded. As described 
in [7], the power multiplier was assembled and its offsets adjusted to zero (thereby setting 
k x &2 and k 3 of (2.1) to zero), such that its output voltage (v a ) would be given by. 

v *=-W x -v y 

where k M (the multiplier gain) is determined from: 

2 .R 

= {from [7] /, » 0.5 mA) 

x y \ 

Substituting /? L - 120 K£2, R - 30 KI2, R y - 62 KI2, and 7,-0.5 mA gives: 

k M = 0.258 

This value of k M will ensure a multiplier output voltage (v fl ) within the output voltage 
range of ±10 V. 


At midband, the feedback resistance (Rj) transforms the final stage into a 
transresistance amplifier, converting it into the current-dependent voltage source 
mentioned in section 2.1.4. This can be deduced from a two port network, midband 
analysis upon the small signal equivalent circuit of the final stage, as shown in Fig. 2.4, 
where again R source represents the internal impedance of the arcjet power supply 



V CONTROLLER 
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Fig. 2.3 Final simulator circuit 
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Z - — - L 

Zl 1 + A-0 


(2.4) 

n 

+b- 

(2.5) 

v c A 

/, 1+A-0 

(2.6) 


where: 


A = - 


R L ~ R D + R SOURCE 

s m ~ |^l 1*0 


P = - 


(2.8) 


(2.7) 


Since R f » R v and usually g m » Rj l , R,» R L then (2.7) can be approximated as: 


As3 ~8m R L R 1 


(2.9) 


If the final stage component values are adjusted such that A-0 » 1, then (2.6) can 
be reduced to: 

(210) 


Substituting for i, gives: 


comp 


, s i 


( 2 . 11 ) 


Thus, by using R f to transform the final stage into a transresistance amplifier, the transfer 
function of the final stage can be made relatively independent of uncontrollable 
parameters which change with output loading (specifically g m and R source)- will 

help to produce a simulator behavior which is independent of the arcjet supply converter’s 
internal impedance, and aid in stabilizing this two-loop feedback system. 


From (2.8) and (2.9), it can be seen that lower values of R f will result in higher 
values of A-0 (thus increasing the validity of (2.11)). But from (2.11) it can be seen that 
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lower values of /^will also reduce the magnitude of the final stage transfer function. This 
can be compensated for by reducing /?, by the same percentage as R f , thereby keeping the 

V 

transfer function magnitude constant. Since i, = — and the current i x is supplied by an 

LM318 op amp, the minimum value of R x will be determined by the output current 
capabilities of the LM318. For this design, /?, was set to 470 Q, the minimum allowable 
value which did not exceed the output current capabilities of the LM318. R f was in turn 
set to 39 k£2, which was high enough to provide sufficient gain for the amplifier, but still 
low enough to ensure that A*(5 » 1. 


R D was chosen to be 4 Q in order to relieve the power dissipation of the final stage 

transistor by approximately 40%, at an operating point corresponding to a simulator 
voltage of 100 V and current of 10 A, which is the expected operating point of an actual 
1 kW arcjet thruster. It was also found that a gate resistance of 120 Q was sufficient to 
prevent parasitic oscillations in the paralled transistors. 


The compensator values were chosen by examining 
simulator. Referring to Fig. 2.1, the gain of the compensator is 

1 r D n 


r * 


= - 


1 + s-R-C a v___ 

a a comp 

S-RyC V 

2 o a 


( 2 . 12 ) 


and, with k x and k^ of (2.1a) equal to zero, the multiplier output volatge is given by: 

v a =/0\ v) = kj-i-v where kj = k M -k.-k v 


However, for the circuit arrangement of Fig. 2.3, the multiplier output voltage (v a ) will be 
negative. Therefore, for the final working version of the simulator, the multiplier output 
voltage is given by: 

v a =/(«', v) = -kfi-v 


Replacing v a by the sum of its DC and incremental parts, and Ayi-v by its Taylor series 



15 


expansion about the operating point, yields: 

'.-v.+o.-zM+if 

op 


• v + H.O.T. 


(2.13) 


op 


Noting that V a =/{/, vj, and assuming small incremental values so that the higher order 
terms (H.O.T.) of (2.13) can be neglected, yields the approximate incremental multiplier 
output voltage: 

v, * - fkr-I-v + kj-V-i ^ (2.14) 

To find the gain of the final stage the transistor must first be replaced with its 
small signal equivalent circuit. Recalling that the transistor actually consists of 16 
paralleled IGBT’s results in the transistor small signal equivalent circuit of Fig. 2.5(a). 
Using PSpice [8] to examine the AC steady state transfer function reveals that the circuit 


16 IGBT 's 




(b) SlnpliFled equivalent circuit 


Ca) Sna t t signal eqlvalent circuit of 
1 6 par a l I ed IGBT 's 


Fig. 2.5 Simplification of transistor small signal equivalent circuit 
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of Fig. 2.5(b) is identical to that of Fig. 2.5(a), as expected. Using the result of Fig. 2.5, 
the small signal equivalent circuit of the final stage can be drawn as shown in Fig. 2.6. 



Fig. 2.6 Final stage small signal equivalent circuit 


The hand analysis of this circuit is somewhat complex. However, as shown in 
Fig. 2.7, a PSpice [8] simulation of Fig. 2.6 shows that it has a dominant-pole frequency 
response for all values of R source ranging from zero to infinity. (For this simulation final 
circuit values of R { R G R f and R D were used while values for C ge C gc C ce and g m were 

estimated from the IGBT data sheet) Therefore, neglecting the effects of higher order 
poles and zeroes, a simplified circuit having a transfer function which very closely 
resembles that of Fig. 2.6 can be obtained. This circuit is shown in Fig. 2.8, where from 
hand analysis it can be shown that the midband gain ( A mUJ ) and dominant high frequency 


pole (c o H ) are given by: 


-8x' R fR L 


Amid tf /+ * f (l Vcomp 


(2.15) 


(O w = 


( /?|1 1 ( T +^ r } C ' (1 +g l’ R L )) 


(2.16) 
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Fig. 2.7(a) Pspice Bode magnitude pi 
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R L “ R D + ^SOURCE 


Fig. 2.8 Simplified circuit having a transfer function closely resembling that of Fig. 2.6 


Usually g x -R L » 1 and (2.15) can be approximated as: 

8i R r R L 


A "*' = 


(2.17) 


and when R L is large enough such that g x R t R x » RA2.ll) reduces to: 


R f 


(2.18) 


which is the same result obtained in (2. 11). 

Since Figs. 2.6 and 2.8 produce similar transfer functions, g x can be found by 


setting 
Therefore: 


A mid 


of (2.15) equal to the A mid of Fig. 2.7(a) for R somCE = 0 (F t = F D ). 


A mid‘( R l +R f) ? ?1 , 


(2.19) 


Similarily, C can be found by setting ( 0 H of (2.16) equal to 2-K-f H of Fig. 2.7(a) for 
r source = °- giving: 

C = = 922.5 pF (2.20) 




Using these values of g, and C along with the final circuit values of R x Rj and R D ’ 
a PSpice [8] frequency response analysis of Fig. 2.8 yields the Bode plots of Fig. 2.9. 
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Comparing Fig. 2.9 with Fig. 2.7 it can indeed be seen that the transfer functions of Figs. 
2.6 and 2.8 closely resemble each other. Therefore, referring to (2.15) and (2.16), the 


gain of the final stage can be closely approximated as: 


Afi( s ) - A m ui’ 


1 


( 2 . 21 ) 


1 +- 


(0 


H 


\ J 

where as given in (2.15) A mid (and therefore A fi (s)) is a negative number. 

The loop gain of the simulator can now be found by breaking the simulator circuit 
at a specific point, as shown in Fig. 2.10, and calculating the ratio of the output voltage to 
the input voltage at this point. 

i 



SOURCE 


Fig. 2.10 Circuit used to calculate the simulator loop gain 


From Fig. 2.10: 


A _ A / 

V c - V comp-^fi' 


,-A fi (s) (2.22) 

" & SOURCE A , A , x ^ SOUR CE 

V V R n + COmP ' A (> Rn + R< 


A A 


v SOURCE 


' SOURCE 


(2.23) 


A 

❖ V 

f 


O' 


comp fs 


Ads)- 


1 SOURCE 


+ R SOURCE 


1 SOURCE 


x SOURCE 


(2.24) 
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^ comp V a 'A. com p(s) 


(2.25) 


•A (s) 

.comp y ' 


Plugging (2.14) into (2.25) gives: 

and substitution of (2,23) and (2.24) into (2.26) gives: 

= - ( kr-IV.AJs).Y — — ) -A (j) 

V comp * *®t*« J ' 


(2.26) 


V = 
comp 


(2.27) 


where, 


y= 


1 SOURCE 


R D + R SOURCE 


(2.28) 


Now, the loop gain can be found from: 


LG — (kfM A U>Y - (2.29) 

V comp V K SOURCE ) 


Substitution of (2.28) into (2.29) and simplifying yields: 

rr a ( \ a , A^t'V'Rsourcf-V)^ 
LG = -A A (s)-A comp (s} 


y R D + R SOURCE j 


(2.30) 


From (2.30) it can be seen that the loop gain becomes positive when I-R source 

Therefore, the following condition must be satisfied for stability: 

(2.31) 


<V. 


R >¥- 

n SOURCE * j 


The result of (2.31) is entirely expected when a load device (such as the simulator) 
having the I-V characteristics of a constant power (i-v = k) curve is connected to the 
output of an arcjet power supply having an internal impedance of R source , as shown in 
Fig. 2.11. 

□ 

^ SOURCE 


V. 


V N/ 



i 

LDAD 





v = f( i ) 


Fig. 2.11 Arcjet power supply having an internal impedance of R source and its load 
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The incremental circuit of Fig. 2.11 is shown in Fig. 2.12 where z represents the 


dynamic impedance of the load which can be found from 

A ^ P 


A 

A V 

Z = ■* = 
l 


dv 

l dl Xv 


dm 


di 


(2.32) 


X.v 


Note from (2.32) that z is also equivalent to the slope of the tangent to the I-V curve at the 
operating point (I,V). 



Fig. 2.12 Incremental circuit of Fig. 2.1 1 


For stability the net resistance of Fig. 2.12 must be positive. Therefore R source 
must be greater than the magnitude of z, or 


R source > I 2 ' 


(2.33) 


Substitution of (2.32) into (2.33) yields: 


R SOURCE > 


df(i) 


di 


Xv 


which for the case of a constant power load 


v =/(i) = 7 


(2.34) 


reduces to (2.31). 


Returning to the simulator loop gain derivation, the substitution of (2.12) and 
(2.21) into (2.30) gives: 


LG = A 


f 1 ' 


mid 


1 +- 


CO 


'H 


1 + s-R-C. 

a a 


s-Rj-C. 

\ 2 a J 


kj.UR 


SOURCE 


-vrt 


R D + ^ SOURCE 


(2.35) 



23 


Therefore, to increase system bandwidth the zero of the compensator was chosen 
such that it would cancel the pole of the final stage. In addition the R 2 C a product was 
made as small as possible to further increase the bandwidth. This product was chosen 
such that the unity-gain crossover frequency of (2.35) would provide maximum 
bandwidth without interference from higher order poles and zeroes of the IGBT’s, op 
amps, and the multiplier. R 2 must also be large enough such that the output current 
capabilities of the LM318 are not exceeded, and small enough such that the LM318’s 
offset currents do not interfere with simulator operation. For this design R 2 was therefore 
chosen to be 22 kO. Then C a was chosen to be 100 pF in order to provide sufficient 
system bandwidth. 

To choose R a , first note from (2.16) that f H can vary from approximately 5 kHz 
(R source = °°) to approximately 30 kHz ( Rsovrce = 0 )- Therefore, to cancel the final stage 
pole, the compensator zero should be placed somewhere in this range. As seen in Fig. 

2.13, the worst case scenario occurs when fa) w < -g In this case, the gain falls off 

A.’C. 

a a 

significantly between <o w and * and the bandwidth is significantly reduced. In 

addition, if the 0 dB crossover should occur on the -40 dB/decade slope, the phase margin 

would be extremely small resulting in a conditionally stable system. To achieve desired 

simulator performance it is important that the simulator be free of unintended internal 

instabilities. Otherwise, possible internal instabilities of the simulator may be falsely 

interpreted as converter instabilities and much time may be wasted before the actual cause 

is determined and corrected. If the simulator is free of internal instabilities there will be 

no doubt as to the cause of any system (the converter and simulator) instabilities. 

Therefore, to minimiz e the scenario of Fig. 2.13, r-=- was placed closer to the lower 

n ‘C_ 

c a 

range of f H (=10 kHz). This results in a value of R a - 150 k£2. 
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Fig. 2.13 Worst case scenario of zero placement 
It should be noted that the reference voltage (v.^) of Fig. 2.3 was supplied by a 

different method than that suggested in the initial circuit design of Fig. 2.1. Both 
methods, however, produce the same large signal, steady state operation equation, given 
by: 



Thus, by adjusting the power potentiometer of Fig. 2.3 (changing v ref ) it is possible to 
shift the I-V curve from one level to another. 

23 SIMULATOR I-V CHARACTERISTICS 

Both the simulator and thruster have bifurcated I-V characteristics. The complete 
simulator I-V characteristic is shown in Fig. 2.14. Initially the final stage transistor will 
be in the off state and no current will flow in the simulator. This corresponds to an 
operating point along I-V charcteristic #1 in Fig. 2.14, and is a simulation of the 
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I 


Fig. 2.14 Complete simulator I-V characteristic 


thruster-off condition. The startup circuit in Fig. 2.3 was added to provide a convenient 
means of starting the simulator. When the startup button is pressed the final stage 
transistor turns on and the operating point jumps from I-V characteristic #1 to a point 
along I-V characteristic #2 in Fig. 2.14, thereby providing a simulation of the thruster-on 
condition. Instabilities of the combined arcjet supply-converter and simulator 
encountered during simulator operation will cause a jump from I-V characteristic #2 back 
to I-V characteristic #1, just as with the actual arcjet thruster. 

For high enough simulator currents the IGBT’s will saturate causing the collector 
to emitter voltage (v ce ) to become nearly zero. As a result, I-V characteristic #2 of Fig. 

2.14 will shift from curve (a) to curve (b), which will have a slope of = R D . Data for high 
current operation of the arcjet is not available. However, it is expected that for high 
currents the I-V characteristic of arc producing devices (such as the arcjet) will have a 
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positive slope, as shown in Fig. 2.15 for the case of TIG welding of aluminum [10]. This 
topic will be discussed further in Chapter 5. For now it is simply noted that the positively 
sloped portion of the I-V characteristic corresponds to current values too high for normal 
arcjet operation. 



Fig. 2.15 I-V characteristics for the TIG welding of aluminum for two different arc 
lengths and gases (Source: [10]) 


2.4 STABILITY REQUIREMENTS FOR VARIOUS ARCJET POWER SUPPLY 
ARRANGEMENTS 

As previously mentioned, the I-V characteristics of the arcjet and of the simulator 
are bifurcated. This bifurcated I-V characteristic is shown again in Fig. 2.16 as the bold 
line curves i-0 and i-v = k. The load lines of an arcjet power supply having an internal 
resistance of R source ^ ^ so shown in Fig. 2.16 for two different values of arcjet supply 

output voltage. As shown for load line #1 of Fig. 2.16, if the arcjet supply output voltage 
is high enough, the load line of the arcjet power supply will intersect the simulator I-V 
characteristic at three points (points (a), (b) and (c) in Fig. 2.16). Point (c) is a stable 
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Fig. 2.16 Simulator I-V characteristic and arcjet supply load lines 

operating point. It corresponds to the condition where the simulator is in its off state, thus 
conducting zero current and having a voltage equivalent to the arcjet power supply open 
circuit output voltage. The stability of points (a) and (b) can be determined by applying 
the stability condition of (2.34) to the respective points in Fig. 2.16. 

From Fig. 2.16, note that the magnitude of the slope of the tangent to the i-v = k 
curve at point (a) (I £1 at point (a)) is greater than the magnitude of the slope of the load 
line ( ^source )• Thus, (2.34) is not satisfied and this intersection point is unstable. 
However, the magnitude of the slope of the tangent to the same curve at point (b) (Izl at 
point (b)) is less than the magnitude of the slope of the load line ( ^source ■)• Thus, (2.34) 
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is satisfied and intersection point (b) is a stable operating point. 

Therefore, for load line #1 of Fig. 2.16, (b) and (c) are both stable operating 
points. Which of these operating points the simulator will reach depends on the dynamics 
which occur after the startup button of Fig. 2.3 is pressed. Repeated testing has shown 
that these dynamics result in simulator operation at point (b) as long as the arcjet power 
supply voltage is sufficiently large such that its load line intersects the constant power 
curve. 


As shown for load line #2 of Fig. 2.16, if the arcjet power supply voltage is not 
sufficiently large then the load line will intersect only with the curve i-0 at point (d). It 
has already been stated that this will be a stable operating point. 

Thus for an arcjet power supply having an internal impedance R S ource » ^ has 
been determined that for stable operation along the I-V curve, the slope of the load line 
must be steeper than the slope of the tangent to the I-V curve, at some point on the curve. 
This is the same as the result stated in (2.34). 

The theoretical discussion of Fig. 2.16 can be demonstrated experimentally by 
using the test set-up shown in Fig. 2.17, where the power supply was adjusted such that 
the load line was close to separating from the power curve, and the peak square wave 
voltage was made sufficiently large to cause the load line to separate from the power 
curve when the square wave reached its negative value. By setting the Tektronix 11401 
Digitizing Oscilloscope up for single trigger mode [9] and closing the switch of Fig. 2.17 
the simulator voltage and current waveforms of Fig. 2.18 were obtained. 

During time interval t x of Fig. 2.18, the load line corresponds to load line #1 of 
Fig. 2.16 and the simulator reaches its stable operating point at point (b). At the end of /, 
the rising square wave voltage causes the arcjet power supply voltage to increase thereby 
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shifting the load line further away from the origin. Thus, the stable operating point will 
move down the I-V curve to the right resulting in a lower simulator voltage and a higher 
simulator current, as shown for time interval t 2 of Fig. 2.18. At the end of t 2 the fallin g 
square wave voltage causes the arcjet power supply voltage to decrease thereby shifting 
the load line to a position corresponding to load line #2 of Fig. 2.16, which is closer to the 
origin. Thus, the operating point moves back up the I-V curve to the left until the load 
line separates from the I-V curve, after which it follows the load line until it reaches its 
final position at point (d) in Fig. 2.11. Note that for time interval t 3 of Fig. 2.18, the 

current is not equal to zero. This is because of the finite time needed to charge the 
compensation capacitor (C a ) to its final value (= 8 ms). 

From Fig. 2.16 it can also be seen that larger values of R^hce will result in more 
stable operating points since transient disturbances of the arcjet power supply would then 
be less likely to separate the load line from the I-V curve. 

For more general cases the internal impedance (Z(s)) of the arcjet power supply 
will include some capacitance and inductance along with the resistance Rso^ce- For this 
more general case the incremental circuit of the arcjet power supply and arcjet is shown in 
Fig. 2.19, where v dist represents a transient disturbance voltage. 



Fig. 2.19 Incremental circuit of the arcjet and arcjet power supply having a more general 
internal impedance of Z(s) 
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The transfer function 

'-v dit "i+Z( S )--r d + Z(s) 

Therefore, if the response to the disturbance voltage (v) is to die out the poles of T must 
be in the left-half plane. Thus, for the more general case the stability requirement can be 
stated as: 

the zeros of {-r d + Z(s)} must be in the left-half plane. 

Since the static I-V curves of the arcjet shown in Fig. 1.2 resemble that of Fig. 
2.16, the above discussion on simulator stability also applies for arcjet operation. 




Chapter 3 

SIMULATOR DYNAMIC IMPEDANCE 

3.1 INTRODUCTION 

If the simulator is to perform as desired, its static and dynamic impedances should 
resemble those of an arcjet thruster as described in Hamley’s paper, "Arcjet Load 
Characteristics" [4]. 

3.2 STATIC IMPEDANCE 

The test set-up of Fig. 3.1 was used to obtain the static impedance data of the 
simulator. This was accomplished by adjusting the power supply voltage in order to 
obtain simulator current and voltage data for simulator currents ranging from 6 A to 15 A, 
at three different settings of v ref . 

As discussed in Chapter 2, R serles of Fig. 3.1 must be made large enough to ensure 
that the power supply load line does not separate from the simulator power curve before 
the static impedance data have been recorded down to 6 A. Also note that for Fig. 3.1: 

V PS = R series ' 1 + V 



Fig. 3.1 Experimental set-up used to obtain static impedance data of the simulator 
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Thus, in order to record static impedance data up to 1-15 A, R series must be small enough 

such that the voltage determined by R series - 15 plus the simulator voltage at 1-15 A and the 

present setting of v ref (V l5 ), does not exceed the output voltage capabilities of the power 

supply (500 V for this case). In other words, R series must also satisfy: 

V w -V, 5 500 -V 15 

R < - = - 

15 15 

By setting R series to 18.3 Q the simulator volt-ampere curves of Fig. 3.2, for several 
settings of v ref , were obtained. For the purpose of comparison, a 1 kW constant power 
curve (shown in dotted lines) is also shown in Fig. 3.2. 

Examination of Fig. 3.2 reveals that the curved portion of the simulator I-V 
characteristic does not actually follow a constant power curve as expected from the design 
discussion of Chapter 2. The only possible explanations for this are that either an error 
was made when zeroing the multiplier offsets (*, and £,) of (2.1) or that some offset 
drift occurred between the time the offsets were adjusted to zero and the time the static 
impedance measurement was made. In either case, the fact that the simulator I-V curves 
are not constant power curves is because the offset coefficients (£, and fcj) of (2.1) are 
not all equal to zero. No steps were taken to correct this because the goal here is not to 
produce simulator I-V characteristics that are constant power curves, but to produce 
simulator I-V characteristics which are similar to those of the arcjet shown in Fig. 1.2. 
Comparing Fig. 1.2 with Fig. 3.2 it can be seen that this goal has been adequately 
accomplished. Also note that the resultant multiplier offset has actually rotated the I-V 
curves of Fig. 3.2 such that they resemble those of an arcjet more closely than those of a 
constant power curve would. 

The effect of changing the power potentiometer setting (changing v re j) of Fig. 2.3, 
is demonstrated in Fig. 3.2 as a shifting of the I-V curve towards or away from the origin. 



SIMULATOR VDLTAGE . 


3 



SIMULATOR CURRENT. A 


Fig. 3.2 Experimental simulator static Volt/ Ampere characteristics for various settings of 
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In addition, the effect of IGBT saturation is also demonstrated as a change from the I-V 
curve to a straight line having a positive slope of * R D . 


3.3 DYNAMIC IMPEDANCE 

In his paper [4], Hamley has shown that the dynamic impedance of an actual arcjet 
thruster has a negative real component and no imaginary component for frequencies 
ranging from 100 Hz up to approximately 50 kHz. Hamley’s results can be seen in Fig. 
3.3 where the arcjet’s upper cutoff frequency (UCF), defined here as the frequency at 
which the phase differs from its low frequency value by 45 degrees, appears to be 
approximately 40-50kHz. 


The theoretical simulator dynamic impedance can be found by applying a current 


test source to the simulator final stage, as shown in Fig. 3.4. From Fig. 3.4 it is clear that 

A 

A V 

the dynamic impedance ( z = -? ) can be found from, 

A A 


Rr 


- I 


(3.1) 


where from Chapter 2, 

A A 


S. = vW'¥) = - (M* + * ».2) 


where. 




(l +s-R-C a ] 

( ] 

sR 2 C a 

V 2 a J 

i+— 


and as in (2.15) and (2.16), 




-8i*r*L 


/? / +/? 1 .(i+* 1 .* i ) 


(3.3) 


<»// = 


Rr 


(*."< T7^y)-(c (1 +*,/y> 

Plugging (3.2) into (3.1) and solving for the dynamic impedance 


/ A \ 
' v ' 

i 

\ J 


, yields: 
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Mass Row [kg/sec] 

O 4.8x10 * 5 
□ 4.5x10 -5 
A 4.3x10 -S 


Fig. 3.3(a) Dynamic impedance magnitude vs. frequency for varoius flow rates (I— 10 A) 



Mass FIow,kg/sec 
O 4.8x10 * 5 
□ 4.5x10 -5 
A 4.3x10 * 5 


5 


Fig. 3.3(b) Dynamic impedance phase angle vs. frequency at various flow rates (I— 10 A) 


Fig. 33 Arcjet dynamic impedance data (Source: [4]) 
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Fig. 3.4 Circuit used to derive the simulator theoretical dynamic impedance 


a v (Ro-k^V-X) 
z = ■* = ■ 


i (1 +k r I-X) 
Substituting for* from (3.3) gives: 


(3.4) 


A 


Ru+kj-V 

(l +sR-C a ) 

a 0 

(A \ 

n mid 

S'RyC 

\ 2 ° ) 

1 + — 
% 

1 - kjb 

l+s-R -C a \ 

rA y 

n mid 

s-Rj'C 

< 2 a ) 

1 +- 2 - 


(3.5) 


If, as discussed in Chapter 2, the compensator zero is made to cancel with the final 
stage pole, then the dynamic impedance reduces to: 

(3.6) 


A R ° + S 'R 2 C c 
Z = - 


1 - 


s R 2 C a 


Rearrangement of (3.6) yields 
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A 

Z = 


*D- 


S + 


^T'^'^mid' 7 


where 7 = 


1 ' 

^C a 


(3.7) 


s ~ MA*r T 

where as shown in (2.15) A mid is a negative number. (3.7) can also be written as: 


A 

Z - 


r 

" / 

\ j 


(0. 


1 + — 
v p J 


(3.8) 


where: 


0 ) = 

z 


KA ^ 7 


Rr 


and “ = A mii kfbi I 


For low frequencies (3.7) and (3.8) reduce to, 

A V 

za ~ 7 


(3.9) 


which as discussed in Chapter 2, is the expected value of dynamic impedance for a load 
device which has a constant power I-V characteristic. Thus, for low frequencies the 
theoretical dynamic impedance of the simulator consist of a negative real component, in 
agreement with that of an actual arcjet thruster. 

Using the HP Dynamic Signal Analyzer (model 3562A) in its synthesis mode 
[11], plots of the imaginary and real parts of (3.7) vs. frequency were obtained for 
simulated load currents of 6, 10 and 14 A. The upper cutoff frequency (UCF) was also 
calculated for each of these operating points and is noted in the upper right comer of each 
plot (note that the measurement range of the HP 3562A is limited to: 64 \iHz - 100kHz). 
The results, shown in Figs. 3.5, 3.6 and 3.7, are consistent with those obtained from an 
actual arcjet thruster as given by Hamley’s paper and shown in Fig. 3.3 [4]. 


Finally, using the experimental setup shown in Fig. 3.8 and operating along static 
I-V characteristic #1 of Fig. 3.2, plots of the imaginary and real parts vs. frequency of the 
actual simulator dynamic impedance for simulator load currents of 6, 10 and 14 A were 
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Fig. 3.5 Real and imaginary components of (3.7) vs. frequency for 1-6 A and V-148 V 
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Fig. 3.7 Real and imaginary components of (3.7) vs. frequency for 1-14 A and V-78 V 
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Fig. 3.8 Experimental set-up used to obtain the dynamic impedance of the simulator 


obtained. In addition, the simulator upper cutoff frequency for each of these operating 
points was measured and noted. These results are shown in Figs. 3.9, 3.10 and 3.1 1. 

As previously noted the low frequency simulator dynamic impedance, at a given 
operating point, is also equivalent to the slope of the tangent to the static I-V curve at that 
operating point. Thus, comparing the low frequency real part values of Figs. 3.9-3.11 
with the slopes of the tangents at points (a), (b), and (c) of Fig. 3.2, shows close 
agreement between corresponding low frequency dynamic impedance values (Note that 
the low frequency dynamic impedance values of Fig. 3.5-3.7 and 3.9-3. 1 1 are given by Y a 

in the upper left hand comer of each plot). Differences can be attributed to the limited 
number of sampling points used to construct the static I-V curve, and to the limited degree 
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Fig. 3.9 Simulator experimental dynamic impedance components vs. frequency for 1-6 A 
and V-148 V 
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Fig. 3.10 Simulator experimental dynamic impedance components vs. frequency for 
1-10 A and V-102 V 
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of accuracy used to construct the tangent at each operating point. 

As shown in Table 3.1, there is a slight difference between the theoretical and 
experimental low frequency dynamic impedance values at each operating point. The 
reason for this difference is that the theoretical low frequency dynamic impedance value 
is given by (3.9), which as already mentioned is equivalent to the slope of the tangent to a 
constant power curve at the operating point (I,V). Whereas, the experimental low 
frequency dynamic impedance value is given by the slope of the tangent to static I-V 
curve #1 of Fig. 3.2, at the operating point (I,V). And as discussed in section 3.2, I-V 
curve #1 of Fig. 3.2 is actually rotated from that of a constant power curve in such a 
manner that it produces a smaller tangential slope magnitude at a given operating point 
than a constant power curve would, which verifies the results shown in Table 3.1. 


! Low Frequency Dynamic Impedance 

Values 

Operating 

1-6 A 

I— 10 A 

1-14 A 

Point 

V-148 V 

V-102 V 

V-78 V 

Theoretical 

-24.67 O 

-10.20 0 

-5.57 O 

Experimental 

-18.56 0 

-8.73 O 

-4.97 O 


Table 3.1 Comparison of theoretical and experimental simulator dynamic impedances 

Overall, the experimental results shown in Figs. 3.9-3.11 are consistent with the 
synthesis analysis (theoretical results), of Figs. 3.5 - 3.7, and Hamley’s results [4]. Thus 
it can be said that the static and dynamic impedances of the simulator closely resemble 
those of an actual arcjet thruster. Note that the simulator upper cutoff frequency 
(=60kHz) agrees closely with that of the arcjet thruster («40-50kHz) for an output current 
of 10A. The synthesis plots show a higher upper cutoff frequency (=100kHz) because the 
synthesis analysis does not take into account any of the simulator’s higher order poles and 


zeroes. 














Chapter 4 


SIMULATOR PERFORMANCE 

4.1 INTRODUCTION 

The purpose of this chapter is to examine simulator performance. This was 
accomplished by comparing the results obtained from testing the simulator with a 
cunrent-mode-controlled (CMC) converter, with a load-current-regulating outer control 
loop designed for 1 kW arcjet applications, with those obtained by testing an arcjet 
thruster with the same converter. The thruster testing was performed at the NASA Lewis 
Research Center using a NASA 1 kW arcjet thruster. For the test set-up, shown in Fig. 
4.1, the CMC converter was adjusted for a constant DC output current of 10 A. The 



Fig. 4.1 Experimental set-up used for testing with a CMC converter 


CMC converter used for this testing is shown in Fig. 4.2. This converter has been 
described in detail elsewhere [12]. The ignition circuit referenced in Fig. 4.2 is shown in 
Fig. 4.3. This ignition circuit was developed following the guidelines set forth in [13]. 
The purpose of the ignition circuit is to provide a high voltage pulse between the anode 
and cathode of the arcjet thruster to ignite the arc. 
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Cout =3000pF 



Fig. 4.3 The ignition circuit referenced In Fig 4.2 
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4.2 STEADY STATE OUTPUT RIPPLE 

The CMC converter of Fig. 4.2 was initially designed to have an output 
inductance ( L ou ) of 384 pH. However, to ensure a wide enough ignition pulse to ignite 

the arc, the converter was also tested with an output inductance of 2 mH. Therefore the 
steady state output ripple results will be given for both values of output inductance. 

4.2.1 Steady State Output Ripple with L 0 ut = 384 \iH 

The steady state output ripple characteristics of an arcjet thruster tested with the 
CMC converter of Fig. 4.1 ( L out = 384 \xH) are shown in Fig. 4.4(a), while those of the 
simulator tested with the same converter are shown in Fig. 4.4(b). The simulator results 
of Fig. 4.4(b) were obtained by adjusting the power potentiometer of the simulator 
(shown in Fig. 2.3) so that the converter DC output voltage was equal to that of Fig. 
4.4(a), which is approximately 90 V. 

Observation of Fig. 4.4(b) reveals that the simulator is indeed performing as 
expected, producing output ripple characteristics which clearly indicate a negative 
simulator dynamic impedance (the change in voltage to the change in current ratio is 
negative for a given time period), as expected from the discussion of Chapter 3. 

Comparison of Figs. 4.4(a) and (b) reveals similarities to the extent that they both 
have a negative dynamic impedance and they both exhibit similar ringing characteristics. 
The causes of ringing will be discussed in section 4.2.2. The figures are different to the 
extent that they appear to have unequal dynamic impedances, and the ringing is more 
damped in the simulator results of Fig. 4.4(b) than it is in the arcjet results of Fig. 4.4(a). 
Calculating the dynamic impedance from the change in voltage and current data shown 
for each figure, yields: 


z = 4r = -3.60 for FigA.Ma) 
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50 V 


div 
2.5 A 


A/®.55A| output voltage 


t output current 


12.5 ps/div 

Fig. 4.4(a) Arcjet steady state output ripple characteristics for Low = 384 ]\H 
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Fig. 4.4(b) Simulator steady state output ripple characteristics for Low = 384 \iH 
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z = ^r~ -6.5ft for Fig A. 4(b) 

Thus, the simulator has a higher value of dynamic impedance than the arcjet. Comparing 
the slopes of the tangents to the I-V curves of Fig. 1.2 with those of Fig. 3.2, for a given 
current, shows that this was an expected result 

The circuit responsible for the output ringing can be simplified to that of Fig. 4.5, 
where L represents the converter stray and/or leakage inductance involved, C om represents 

the ignition capacitor of the converter, r d represents the dynamic impedance of the 
simulator, and R hss represents the total losses of the resonant circuit of Fig. 4.5. 


L 



Fig. 4.5 Simplified output ringing circuit 


If < r d (usually the case for Fig. 4.5) and R eq is defined as the parallel 


combination of r, and that is. 


n _ r d R u™ 
— — 

49 ratios. 


and remembering that r d is a negative number, it becomes apparent that higher 
magnitudes of r d will give lower values of R t resulting in more damping of the output 
ringing. Therefore, the fact that Figs. 4.4(a) and (b) have unequal damping can be 
attributed to their unequal dynamic impedances, and as expected the figure with the 
higher magnitude of r d (Fig. 4.4(b)) has the greater damping. 
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Later (in Chapter 5), it will be shown how to change the dynamic impedance of 
the simulator so that it more closely resembles that of the arcjet thruster. For now, it is 
simply noted that the simulator and thruster produce similar output ripple characteristics 
when tested with a CMC converter having an output inductance of 384 \iH. 

4.2.2In vestigation of Output Ringing 

A plot of the CMC converter transformer primary current (when tested with the 
simulator) is shown in Fig. 4.6. During the time interval t a shown in Fig. 4.6, the 
appropriate pair of bridge transistors (Q1 & Q4 or Q2 & Q3) along with the appropriate 
rectifier diode (D1 or D2) of Fig. 4.2, are conducting. Therefore, for the time interval t a , 
the incremental converter circuit of Fig. 4.7 is in effect 

Examination of Fig. 4.7 reveals that the ringing during time interval t a , of Fig. 4.6, 
is caused mainly by the primary to secondary leakage inductance of the transformer, 
referred to the secondary, acting in conjunction with the 3000 pF ignition capacitor (C out ). 

A plot of the CMC converter transformer secondary current (when tested with the 
simulator) is shown in Fig. 4.8. During the time interval t b shown in Fig. 4.8, the 
appropriate pair of bridge transistors (Q1 & Q3 or Q2 & Q4) and both rectifier diodes (D1 
& D2) of Fig. 4.2, are conducting. Therefore, for the time interval t b , the incremental 
converter circuit of Fig. 4.9 is in effect 

Examination of Fig. 4.9 reveals that both the transformer primary and secondary 
windings are directly shorted. Therefore the ringing during time interval t b , of Fig. 4.8, is 
caused mainly by the secondary to secondary leakage inductance plus the stray inductance 
of the secondary circuit, acting in conjunction with the 3000 pF ignition capacitor ( C out ). 
This conclusion is supported by the fact that there is no ringing in the transformer primary 



Fig. 4.7 Incremental converter circuit for time interval t t 
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current of Fig. 4.6 during the time interval t b . 

The ringing during time interval t b has a higher frequency than that of time 
interval t a . This comes as no surprise since the secondary to secondary leakage inductance 
plus the stray inductance of the secondary circuit is expected to be much less than the 
primary to secondary leakage inductance of the transformer referred to the secondary, and 
the ringing frequency is determined from: 

“ = 7Fc 

4.2.3 Steady State Output Ripple with L 0 ut = 2 mH 

The steady state output ripple characteristics of an arcjet thruster tested with the 
CMC converter of Fig. 4.2 ( L out = 2 mH) are shown in Fig. 4.10(a), while those of the 

simulator tested with the same converter are shown in Fig. 4.10(b). In this case, the 
simulator results were obtained by setting the power potentiometer of the simulator to 
obtain a converter DC output voltage equal to that of Fig. 4.10(a), which is approximately 
99 V. 

As expected, the results of Figs. 4.10(a) and (b) are very similar to those of Figs. 
4.4(a) and (b). The only difference between the two sets of figures is a smaller ripple 
amplitude for Figs.4. 10(a) and (b). This is to be expected since the converter output 
inductance ( L out ) is significantly larger for these two figures than it is for Figs. 4.4(a) and 

(b). 

All other aspects of Figs. 4. 10(a) and (b) are the same as discussed in section 4.2. 1 
with Figs. 4.4(a) and (b) replaced everywhere by Figs. 4.10(a) and (b), respectively. 
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Fig. 4.10(a) Arc jet steady state output ripple characteristics for Lout = 2 mH 
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Fig. 4.10(b) Simulator steady state output ripple characteristics for Lous = 2 mtf 
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4.5 LARGE SIGNAL TRANSIENT RESPONSE WITH L ou t =2 mH 

The large signal transient response of an arcjet thruster tested with the CMC 
converter ( L om = 2 mH) is shown in Fig. 4.11(a), while that of the simulator tested with 
the same converter is shown in Fig. 4.1 1(b). These measurements were taken by injecting 
a pulse waveform into the control circuit of the converter, as shown in Fig. 4.12. The 
pulse magnitude was set to the maximum allowable value which did not cause simulator 
or thruster turn-off. 

It is apparent that the current and voltage waveforms of the arcjet shown in Fig. 
4.11(a), closely resemble those of the simulator shown in Fig. 4.11(b). The slight 
differences are believed to be caused by the unequal dynamic impedances of the simulator 
and arcjet. Finally, as expected for a negative dynamic impedance, the current pulse goes 
negative for a corresponding positive pulse in voltage. 

Further examination of Figs. 4.11(a) and (b) reveals evidence of underdamped 
ringing on the current and voltage pulse waveforms. This underdamped ringing indicates 
the presence of a complex conjugate pole pair in the system (the converter and its load) 
transfer function, which is located in the left half s-plane, fairly close to the jco-axis. 
Therefore, the bandwidth and phase margin of the system are relatively small, resulting in 
a marginally stable system. This occurs because the converter loop gain was originally 
compensated by using a resistive load bank in place of the simulator (before the simulator 
was designed). This is demonstrated in Fig. 4.13, where the large signal transient 
response of a resistive load bank tested with the CMC converter, for the same conditions 
as Figs. 4.11(a) and (b), is shown. Examination of this figure reveals that the system is 
more than adequately compensated for a resistive load. 

Thus, to increase system stability, the converter loop gain should be compensated 
when the converter is operating with the actual load it was designed for. Hence, another 
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Fig. 4. 1 1 (a) Arcjet large signal transient response for Lout - 2 mH 
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Fig. 4.13 Resistive load large signal transient response for L out = 2 mH 

useful application of the simulator arises in compensating the converter for use with an 
actual arcjet thruster. 

4.6 LARGE SIGNAL TRANSIENT RESPONSE WITH L 0 «/ = 384 \iH 

The large signal transient response of the simulator tested with the CMC converter 
(L ou t = 384 \xH) is shown in Fig. 4.14. Unfortunately, no corresponding data was taken 
with the arcjet thruster. 

A comparison of Figs. 4.14 and 4.11(b) reveals that the system is much more 
stable with an output inductance of 384 pH than one of 2 mH. It is known that the CMC 
converter loop gain is dependent upon the value of output inductance used [14], 
Therefore, the converter should be compensated for the final value of output inductance 
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Fig. 4.14 Simulator large signal transient response for L out = 384 jj H 


used. 

With the exception of more ripple in Fig. 4.14 than in Fig. 4.1 1(b), the two figures 
are similar in all other aspects. 

4.7 SIMULATOR EXTINCTION TRANSIENT 

The simulator extinction transient was obtained by increasing the magnitude of the 
pulse waveform injected into the control circuit of Fig. 4.12 until simulator extinction was 
achieved. Using the Tektronix 11401 Digitizing Oscilloscope in the single trigger mode 
[9], the output voltage and current waveforms of Fig. 4.15 were recorded. 

Fig. 4.15 shows that after the pulse occurs the simulator current and voltage 
approach a new equilibrium point However, the new equilibrium point turns out to be 
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100*3 'd i_v 

Fig. 4.15 Output voltage and current during simulator extinction 


unstable and the simulator’s current goes to zero and its voltage to the CMC converters 
open circuit value. Thus, the operating point jumps from I-V characteristic #2 to I-V 
characteristic #1 of Fig. 2.14. The large spike in the simulator voltage is caused by the 
rapid change in the current through L M as the simulator current goes to zero. 

Overall, the simulator extinction response seems to accurately predict the arc 
extinction response of an actual arc jet thruster. 
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Chapter 5 

SUMMARY FOR PART I 
5.1 DYNAMIC IMPEDANCE ADJUSTMENTS 

As suggested in Chapter 4, another issue to be considered is how to make the 
dynamic impedance of the simulator adjustable, so that it can be made to more closely 
resemble that of an arcjet thruster. As mentioned in Chapter 2, the dynamic 
impedance of the simulator is equal to the slope of the tangent to the static I-V 
characteristic at a given operating point. Thus, by adjusting the static I-V curves the 
simulator dynamic impedance can be changed. 

A simple way of doing this is by connecting a voltage feedback network to 
the positive terminal of the PI compensator op amp shown in Fig. 2.3. This network is 
shown in Fig. 5.1 where previously existing components are shown in dotted lines, 
and R1 and R5 represent the power potentiometer of Fig. 2.3. 

to nut t ipl itr 

t 

♦ v » 

R2*22kn 



Fig. 5.1 Voltage feedback network used to adjust the simulator dynamic impedance 



65 


As stated in Chapter 2, for large signal, steady state operation, the compensation 
capacitor (C a ) of Fig. 5.1 will block current flow through its path. Thus, the large signal. 


steady state operation of Fig. 5.1 will be determined by, 



where Vj is given by: 


V *4 

V ‘ V /? po( + 200 - 10 3 + /? 4 


Rearranging (5.1) yields: 



(5.2) 


(5.3) 


Recalling from Chapter 2 that v a = - kfi-v (assuming that k x and of (2.1) are all zero) 
and that /? 3 = R 2 (see Fig. 2.3), (5.3) can be reduced to: 

V re/~ 2' v l + kfi’ V (5.4) 


Substituting (5.2) into (5.4) yields: 

2/? 4 


Vref R + 200- 10 3 + * 4 V + fcr < V 


Thus, the constant power I-V curves produced by (2.2a) are being altered by the addition 
of a term proportional to the simulator voltage (v), hence the term voltage feedback 
networic for the circuit of Fig. 5.1. 

By using the existing value of R3 (22 k£2), setting R4 to 4.7 k£2, and adjusting /? pot 
to 1 M£2, in Fig. 5.1, static I-V characteristic (c), shown in Fig. 5.2 (where VF - Voltage 
Feedback), was obtained. For comparison, the static I-V characteristic of an arcjet (the 
middle curve of Fig. 1.2) and the previous I-V characteristic of the simulator without 
voltage feedback (curve #1 of Fig. 3.2) are shown in Fig. 5.2 as curves (a) and (b), 
respectively. Using the values of R3 and R4 mentioned above, and reducing R ^ yields 
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Fig. 52 Static I-V characteristics obtained using the voltage feedback network of Fig. 5.1 
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IV curve (d) of Fig. 5.2. Reducing R to zero yields I-V curve (e). 

Comparing Fig. 5.2 curves (a), (b) and (e) reveals that indeed the dynamic 
impedance of the simulator is being adjusted closer to that of the arcjet. However, curve 
(e) needs to be shifted away from the origin in order to achieve similar operating points as 
the arcjet curve (a). As discussed in Chapter 2 and demonstrated in Chapter 3, this can be 
accomplished by adjusting the power potentiometer of Fig. 2.3. Adjusting this 
potentiometer for an operating point of 1 16.2 V at 10 A (the approximate arcjet voltage of 
curve (a) at 10 A) results in I-V curve (f) of Fig. 5.2. 

Comparing curves (a), (b) and (f) of Fig. 5.2 reveals that the simulator dynamic 
impedance has been decreased, but is still greater than that of the arcjet. These reductions 
have occured as the positive terminal voltage of the LM318, of Fig. 5.1, has been steadily 
increased by reducing towards zero. Thus, further dynamic impedance reductions 
can be made by increasing R4 and again adjusting towards zero. Changing R4 to 10 
k£l, adjusting to zero and adjusting the power potentiometer (R^^) to its maximum 
value (such that R 1-0 in Fig. 5.1) yields I-V curve (g) of Fig. 5.2. Examination of this 
curve reveals that it has a dynamic impedance very close to that of the arcjet curve (a). 
However, curve (g) must be shifted away from the origin to achieve similar operating 
points as (a), and the power potentiometer (R^*.) is already set to its maximum value. 
However, by changing R3 to 10 kQ, the range of R^^ can be extended. Thus, by 
adjusting R^^ for an operating point of 116.8 V at 10 A, curve (h) of Fig. 5.2 was 
obtained. 

Comparing curves (a) and (h) of Fig. 5.2 reveals close agreement between the I-V 
characteristics (and therefore the dynamic impedance) of the arcjet and simulator, 
especially in the range of 8 - 10.5 A. 
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Note that by letting, 


*3 = 


2R. 


+ 200-10 +/? 4 


(5.5) can be rearranged to yield: 


V 

v = L - 


i tj + 1 kj-i 


(5.6) 


which is the same as (2.1b) (repeated here for convenience) with ifc, = k^ = 0: 

(5.7) 


V ref- k 1-M 

V=f^l) = — L 


k} + kj-i 


Thus, the same results achieved by the voltage feedback network of Fig. 5.1 could 
also be achieved by manipulating the multiplier offset potentiometers so that k^ of (2.1b) 

is varied while k x = k^ = 0. This could be accomplished by varying the I offset pot of Fig. 
2.3. 


If in addition to k^, k x and k^ are also varied (by varying the P and V offset pots of 

Fig. 2.3), the degree of freedom available in matching the simulator I-V characteristics to 
those of the arcjet will increase. However, since this will involve the adjustment of three 
potentiometers it is suggested that a curve tracer be employed to give a continuous trace 
of the simulator static I-V characteristic as each pot is varied, so that the best 
approximation to the desired I-V curve can be rapidly found. 

5.2 OTHER APPLICATIONS 

One application of interest is in the area of welding power supply testing. As 
illustrated in Fig. 2.15, the I-V characteristic of a welding arc is similar to that of the 
simulator shown in Fig. 2.14 . 

Using the techniques of the previous section along with a change in R D of Fig. 2.3, 
the simulator static I-V curve can be made to closely resemble that of a given welding arc. 
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In this way the simulator could be used as a load when testing welding power supplies . 

Another possible application of interest concerns the simulation of switching 
converters in solar power systems testing (such as the case discussed in [IS]). When a 
switching converter is in regulation the current will decrease (increase) for a 
corresponding increase (decrease) in voltage such that its output power remains constant 
Thus, while in regulation the switching converter follows a constant power I-V 
characteristic. When the switching converter voltage gets sufficiently low the converter 
will drop out of regulation and behave as a resistive load. Thus, as given in [15], the I-V 
characteristics of a regulated switching converter are as shown in Fig. S.3. 



Fig. 5.3 I-V characteristics of a regulated switching converter (load lines 1-5) 
(Source: [15) 


It has already been shown that the simulator can reproduce the constant power 
portion of the I-V characteristic shown in Fig. 5.3. By changing the sign of the multiplier 
(switching the current sense leads on the emitter resistor) and of the DC voltage applied to 
the power potentiometer of Fig. 2.3, the straight line portion of the simulator I-V 
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characteristic shown in Fig. 2.14 can be made to resemble that of the switching converter 
shown in Fig. 5.3 

By using the techniques of section 5.1 and the appropriate value of R D , the 
simulator I-V characteristic can be made to closely resemble the input characteristic of a 
switching converter having a constant load. Thus, in this way the simulator could be used 
to simulate actual switching converters when testing solar power systems. 

In short, the simulator has possible applications in any situation where the load 
current increases (decreases) with decreasing (increasing) voltage. 

S3 SIMULATOR LIMITATIONS 

As shown in Figs. 3.9-3.11 the simulator upper cutoff frequency is limited to 
approximately 40-60 kHz depending on operating point (Remember that the upper cutoff 
frequency has been defined in this study as the frequency at which the phase of the 
dynamic impedance varies from its low frequency value by 45 degrees.) The limited 
arcjet data obtained from [4] suggests that the arcjet also has an upper cutoff frequency of 
approximately 40-60 kHz. Therefore, it appears that the upper cutoff frequency of the 
simulator is comparable to that of the arcjet. However more arcjet data should be 
obtained before this result is verified. For most applications the simulator upper cutoff 
frequency should be large enough to include the source ripple frequency, which for this 
case is 40 kHz. 

Due to a lack of mathematical equations describing arcjet behavior and a lack of 
sufficient arcjet experimental data to compare with simulator results, the extent to which 
the simulator matches the large signal dynamics of an actual arcjet remains relatively 
unknown. 
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The ignition pulse created by the ignition circuit of Fig. 4.3 will be in the range of 
1-2 kV. Therefore, the ignition circuit and simulator were not tested together as such 
testing would require very high voltage simulator devices. Thus, the simulator should not 
be used for ignition circuit testing. 

5.4 CONCLUSION 

The simulator as developed herein has been shown to produce I-V characteristics 
which closely resemble those of an arcjet thruster, thereby demonstrating its ability to be 
used as a load device in determining the stability of converters designed for arcjet 
applications. In addition, the simulator has been shown to be a fairly accurate predictor of 
actual arcjet performance. Since the system (the converter and its load) loop gain has 
been shown to be dependent on the type of load used, another useful application of the 
simulator arises in compensating the converter for use with an actual arcjet thruster. It 
has also been shown that the simulator can be constructed in the laboratory using 
relatively inexpensive power electronic components. 

The stability of various arcjet supply arrangements has been demonstrated by 
proving that the internal resistance of the arcjet supply converter must be greater in 
magnitude than the dynamic resistance of the simulator to achieve a stable, non-zero 
current, operating point 

The static and dynamic impedances of the simulator were shown to agree closely 
with results obtained from an actual arcjet thruster. In addition, simulator performance 
was determined by testing the simulator with a current mode controlled (CMC) converter 
and comparing the results with those of an arcjet tested on the same converter. The data 
obtained showed close agreement between simulator and thruster performance. 


Methods of varying the simulator dynamic impedance so that it would more 
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closely resemble that of an arcjet thruster, or some other device, have also been discussed. 
Finally, it has been shown that the simulator may have applications in other areas, such as 
in testing power supplies for welding or for the simulation of switching converters in 
testing solar power systems. 


/ 



73 


PART II 

ARCJET IGNITION CIRCUIT DEVELOPMENT 

Chapter 6 
INTRODUCTION 

Fig. 6.1 shows the lkW converter power circuit which was developed in [12]. 
The gate drivers, gate drive current booster, current mode controller (CMC), diagonal 
to phase-shift conversion logic, and CMC Hall effect current sensor are shown in 
Figs. 6. 2-6. 6, respectively. All of this circuitry was developed as part of the study 
described in [12], and the circuit details are described in that reference. The ignition 
circuit is shown in Fig. 6.7, and its theory of operation is described in the next 
chapter. 
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Gate drive circuit 


Fig. 6.2. Isolated gate drive circuit with power-fall cutoff. 
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Fig. 6.4 Current Mode Controller 
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Fig. 6.5 Logic for converting 
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Fig. 6.6 CMC sensor using a Hall-device dc current sensor. 
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Chapter 7 

THEORY OF OPERATION FOR IGNITION CIRCUIT 

This arcjet ignition circuit is shown in Fig. 6.7. This circuit is designed to 
operate from the primary input voltage to the arcjet converter, in the range of 20-35 
V, and provide a consistent ignition pulse to the arcjet. The converter output 

inductor is provided with a closely-coupled extra winding through which energy is 
introduced to this inductor to form the ignition pulse. When an ignition command is 
received by the control portion of the circuit, the IGBT is turned on, allowing the 

current in the extra inductor winding to ramp up to a predetermined set point, at 
which time the IGBT is turned off. This begins an interval in which the inductor 

primary voltage rises resonantly to the breakdown voltage of the arcjet. Assuming 

enough energy is initially stored in the inductor, an initial current is established in 
the arcjet corresponding to an arc voltage below the open-circuit voltage of the 
arcjet converter. This permits the current in the inductor to rise to the point at 

which the converter current-regulation loop takes control of the converter. The 
circuit of Fig. 6.7 is similar to that of [13] in most respects, and was derived from it. 

The control portion of the ignition circuit uses a TSC 429 gate-driver to gate 
Ql, the IGBT. This gate driver is made self-latching by the addition of Q2. An "off' 
state for this latch is assured upon power-up by Q4 and its associated circuitry. An 

ignition pulse is manually triggered by the NO/NC push button shown, together with 
Q3. An externally-commanded ignition pulse could be coupled in through an opto- 
isolator at this point. 

An ignition command sets Ql's latching gate driver, allowing the current in 
the extra winding on inductor L ou t to ramp up at a rate determined by the primary 
input voltage. This ramping current is sensed by the current transformer and Q5, 
and when it reaches its set point the gate driver is reset by Q5. Because a current set 
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point is used instead of a fixed inductor-charging time, a consistent amount o 

energy is stored in the inductor regardless of varying supply voltage. Upon tun 

off of Ql, the inductor voltage rises sinusoidally, the resonance being determined b] 

the 2.22 pF capacitor in parallel with the IGBT and the 3000 pF on the inductor mail 

winding. When break over of the arcjet occurs at the angular time "0 M at the natura 

resonant frequency, an initial current equal to the charging-current set point time 
cos 0 , reflected through the inductor turns ratio, is established in the arcjet. (It i 

generally found that cos 0 is close to one for practical purposes in this design.) It i 

then expected that the arcjet supply converter will establish an appropriati 
operating point for the arcjet before the energy stored in the inductor is dissipated. 

The open-circuit voltage produced by the ignition circuit is equal to its se 
point current multiplied by the characteristic resistance of the resonant circuil 

approximately 3.4 kV with a 2-mH output inductor. Should ignition fail to occur, th 

diode in anti-parallel with Ql allows the return of the stored energy to the 1000 pi 

supply capacitor, and reset of the ignition circuit. 

The method used here is to initially ignite the arc with a high voltage pulsi 

and then sustain the steady state arc with a current source. This process can bi 
explained as follows. The simplified power circuit is shown in Fig. 7.1. S is initiall; 

closed, and ii ramps up to some value 1 1 » at which time S is opened. The reflectei 

value of Ii is Ip, which is the initial value of iL. The equivalent circuit for this perioi 
is shown in Fig. 7.2, and waveforms for u 0 and iL with no arc are shown in Fig. 7.3. A 
u 0 increases, the voltage across the arc gap travels along the vertical axis from O to / 
in Fig. 7.4 until ignition occurs at v 0 = V\. As the arc is established, i D begins to flow 
and the operating point moves to the v 0 vs. i 0 curve for the arc at some point B. N< 
attempt is made to model the exact behavior of v Q vs. io during ignition, but v 0 wil 
quickly decrease as the energy stored in L 0 is dissipated by the arc. Thus thi 
operating point will move down the v 0 vs. i Q curve as vo decreases. The curren 
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regulated converter now has a 100% duty cycle since i 0 has not yet reached the I 0 reg 
regulation value. After the operating point passes point C, the converter will remain 
at 100% duty because i 0 < I 0 re g As io continues to increase beyond point C, v 0 will 
continue to decrease as the output operating point moves down the arc curve. The 
voltage difference, v s - v 0 , appears across L 0 . I 0 reg represents the regulated current 
value that is controlled by the current feedback loop. Because of the delay caused by 
the response time of this loop, io proceeds past Io reg at point D and on to Io max 
which is the absolute open loop limit fixed by the controller. When io reaches Io max 
at point E, the converter begins to operate in the PWM mode to limit i 0 = Io max- As the 
current control loop begins to respond, the duty cycle will decrease to reduce i 0 = I 0 
reg, and the steady state operating point will move back to D. 

The simplified circuit in Fig. 7.5 provides an explanation as to why a current 
regulated supply should be used instead of a voltage regulated supply. As seen from 
Fig. 7.4, at least over part of its operating range, the arc can be approximately 
modelled by a nonlinear R(i 0 ) and va(i G ) in series. Therefore vo will drop as io 
increases. If voltage regulation is used, the equivalent source resistance, Rj, is quite 
small and Ri-Ra(io) < 0, which is an unstable condition. If current regulation is used, 
R1 is very large and Ri - R a (io) > 0. which provides stability. 









Fig. 7.2. Equivalent circuit after 
□pen 0 t=0 and before ignition. 










NAS?_3.DWG 

f or 

0<t<ti 








Lo 



1 

Vo 

- v s 

+ ^ P 


sin wt 

, w- , 



\ 

Co 



o 

C_J 

0 

_J 

1l = 

! P 

COS ' 

rft 




Vi . 

Ii 

- VO 

ltage and 

current at ignition 

F 

iP • 

7.3. 

Vo 

and 

1l 

w i thoui arc . 




■ignition path 


vo vs. io curve f on 
cunnent negulated power s 


la N 


omox 


NAS7_4 .DWG 


7.4. S implified charac t er i st i c curves 
F or . arcjet and power supply. 




Reg u I a t or Arc 


NAS7_5 . DWG 

Vo = Va ( i o ) - i O' Ra ( i o ) 

V I = i o ( R l - i o ) ) f V a ( i o ) 

Ra ( i o ) and Va (io) are nonlinear 


Fig. 7.5. Simplified equivalent circuit for 
regulator and arc. 





Chapter 8 

EXPERIMENTAL RESULTS 


To verify the operation of the converter and ignition circuits, a series of 
tests were performed at the NASA Lewis Research Center where the power supply was 
operated with an arcjet load. 

Fig. 8.1 shows the open circuit ignition pulse without the 100 V. output and 
with Lo = 2mH. Fig. 8.2 shows the pulse under the same conditions but with the arcjet. 
An expanded view is shown in Fig. 8.3. 

Fig. 8.4 shows v 0 and i 0 at ignition with the peak current limit set at about 
40 A. Fig, 8.5 shows these same waveforms but with an expanded voltage scale. This 
peak value of 40A. was judged to be excessive and raised some concern about damage 
to the arcjet nozzle. Therefore, the peak current limiter setting was lowered, and 
Figs. 8.6 and 8.7 show the results for peak limits of about 12. 5A. and 21 A., 
respectively. Fig. 8.8 shows the steady state ripple in v G and i c , while Fig. 8.9 shows 
the response of v Q and i Q to a step change in the current reference voltage. 

In an effort to save weight, it would be advantageous to reduce the size of the 
output filter choke, but there is some concern that ignition may not occur if the 
choke is too small. This concern is based on the fact that a smaller choke will sustain 
the ignition pulse for a shorter period of time. To investigate the behavior with a 
smaller choke, a series of tests were performed using a 384 pH choke in place of the 
original value of 2mH. The resulting waveforms are shown in Figs. 8.10-8.12. Fig. 
8.10 shows that the ignition pulse is now in excess of 4kV, and ignition occurs in 
about 1.25 ps. Fig. 8.2 indicates that with the 2 raH choke ignition occurred at about 
1.5 kV and within about 3.75 ps. Both chokes appeared to be satisfactory since the arc 
was ignited numerous times with each choke, and no ignition failures were observed. 
Fig. 8.11 shows the output voltage and current ripple. The output voltage and current 
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after ignition with I limit : 25 A. peak are shown in Fig. 8.12. The peak values of v 0 
and i 0 in Fig. 8.12 are considerably lower than those in Fig. 8.10 since the limited 
sampling rate of the digital storage oscilloscope apparently missed the peak values 
when used on the longer time scale. 
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Fig. 8.1. Ignition pulse without arcjet and lOOVdc turned off (Lo 
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Fig. 8.3. Expanded view of Fig. 8 - 2 , 15.4us after initial pulse 
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Fig. 8.4. vo and io at ignition with A. peak 
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Fig. 8.5. Same as Fig. 8.4 but with expanded voltage scale. 
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Fig. 8.6. vo and io at ignition with I,. . — 12.5A. peak. 
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Fig. 8.7. vo and io at ignition with Iii ro it~ 21 A. peak 
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Fiq. 8.8. Steady state vo and io (output voltage and current ripple) 
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Fig. 8.9. Response of vo and io to a step change in the current reference 
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Fig. 8.11. Steady state vo and io with 384uH choke (output voltage and current ripple). 
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Chapter 9 

SUMMARY FOR PART II 

It appears that the arcjet power supply functions properly and no operating 
problems were observed. Every attempted ignition was successful, and the control 
circuit appears to be immune from EMI and instability problems. Of course, further 
optimization is always possible, and there are certain areas of investigation that 
might lead to reductions in cost and/or weight. One of those may be a further 
reduction in the output filter choke, since values below 384 pH were not tested. 
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